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Experimental searches for lepton-flavor- violating coherent muon-to-electron conversion in the field 
of a nucleus, have been proposed to reach the unprecedented sensitivity of 10"^" — 10~^* per stopped 
muon. At that level, they probe new interactions at effective-mass scales well beyond 1000 TeV. How- 
ever, they must contend with background from ordinary bound muon decay. To better understand 
. . . the background-spectrum shape and rate, we have carried out a detailed analysis of Coulombic- 

' bound-state muon decay, including nuclear recoil. Implications for future experiments are briefly 

, discussed. 
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I. INTRODUCTION AND MOTIVATION 

From the observation of neutrino oscillations, we now know that lepton flavors (electron, muon and tau number) 
are not conserved. However, the mixing and small neutrino mass differences seen in oscillations have negligible 
effect on charged-lepton flavor violating (CLFV) reactions such as /.t — 67 or r — )■ (now predicted to occur with 



^p^i unobservable tiny branching ratios of about 10 ^"'). So, if charged-lepton flavor violation were to be experimentally 

f*) 1' detected, it would have to come from "new physics" such as supersymmetry, heavy neutrino mixing, leptoquark 
interactions or some other extension of the Standard Model. In that way, charged-lepton number violating reactions 

1 ' provide a discovery window to interactions, beyond Standard Model expectations, reaching effective-mass scales above 

; 0(1000 TeV) |lj,i2|]. 

O^l ' Because muons can be copiously produced at accelerators and are relatively long lived (2.2/is), they have been at 

' the forefront of searches for CLFV [J Q ■ One reaction that can be probed with particularly high sensitivity is the 

' muon-electron conversion in a muonic atom, 

in 



^i- + {A,Z)^e- + {A,Z), (1) 



' where {A, Z) represents a nucleus of atomic number Z and mass number A. Various experiments have been performed 
. over the years to search for this process [p. The most recent, and stringent, results come from the SINDRUM II 
I ' Collaboration Q, which reports an upper limit of 7 x 10^^'^ for the branching ratio of the conversion process relative 
• • , to muon capture in gold, and a similar unpublished bound for titanium. Several new efforts are being planned. In 
^ ' the nearest future, the DeeMe Collaboration [5| has proposed to reach 10"'^'* sensitivity. Larger scale searches, Mu2e 
'kA : at Fermilab [| and COMET at J-PARC 0, aim for sensitivities below 10 In the long run, intensity upgrades 
5_j ■ at Fermilab and the proposal PRISM/PRIME at J-PARC may allow them to reach 10~^^ sensitivity, a limit only 
\ accessible with muon-electron conversion in nuclei. For comparison, the current best upper bound on the branching 
ratio of the decay /i — > 67, set by the MEGA experiment, is 1.2 x 10^^^ (90% confidence level) Q. For some 
mechanisms of CLFV, the conversion process is less sensitive than /i — >■ 67 by a factor on the order of a few hundred 
Q. But even in those cases, a 10^^^ conversion search is more sensitive than this best current bound, and may be 
competitive with the new search for /i — > 67 by the MEG experiment [lo| . In addition, the conversion process is also 
sensitive to CLFV chiral-conserving amplitudes that do not contribute to /x — )■ 67. 

The success of the conversion searches depends critically on control of the background events. The signal for the 
/i — e conversion process in Eq. ([IJ is a mono-energetic electron with energy -E^e, given by 

E^e = — Eh — Ei-cc, (2) 

where is the muon mass, -Eb — 2^^a^m^/2 is the binding energy of the muonic atom, and E'rcc — 'Ti^/(2mAr) is 
the nuclear-recoil energy, with a the fine-structure constant and ttin the nucleus mass. The main physics background 
for this signal comes from the so-called muon decay in orbit (DIO), a process in which the muon decays in the normal 
way, i.e. fi~ e~Vey^, while in the orbit of the atom. Whereas in a free muon decay, in order to conserve energy 
and three-momentum, the maximum electron energy is m^/2, for decay in orbit the presence of an additional particle 
(the nucleus), which can absorb three- momentum, causes the maximum electron energy to be -E^e- Therefore, the 
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high-energy tail of the electron spectrum in muon decay in orbit constitutes a background for conversion searches. A 
detailed study of that background is the main focus of this work. 

Several theoretical studies of the muon decay in orbit have been performed, starting with Ref. 11] about 60 years 
ago. Reference fl^l presented expressions which allow for a calculation of the electron spectrum including relativistic 
effects in the muon wavefunction, the Coulomb interaction between the electron and the nucleus and a finite nuclear 
size. Nuclear-recoil effects which, as we will discuss later, need to be considered in the high-energy region were only 
included in the Born approximation (that is, using the non-relativistic Schrodinger wave function for the muon and 
a plane wave for the electron), which is not adequate for the high-energy tail. Later, Refs. fl3, 14] presented similar 
expressions for the electron spectrum completely neglecting nuclear-recoil effects, evaluating it for several different 
elements. None of these references focused on the high-energy endpoint of the spectrum, which is the region of interest 
for conversion experiments. References [H, HBl did study the high-energy end of the electron spectrum, and presented 
approximate results which allow for a quick rough estimate of the muon decay in orbit contribution to the background 
in conversion experiments. However, a detailed evaluation of the high-energy region of the electron spectrum is still 
missing in the literature. What is typically done, to account for the background from muon decay in orbit, is to 
connect (in a somewhat arbitrary way) the approximate expressions given in Ref. fld\ with the numerical results 
presented in Ref. 14]. Since this is the main source of background @, Q for the oncoming conversion experiments, a 
more detailed analysis is highly desirable. In this work we discuss all the relevant effects that need to be included in 
the high-energy region of the spectrum and present a precise evaluation of it. Our results for an aluminum (Z = 13) 
nucleus (the intended target in Mu2e and COMET) are presented in Figure [T] 

The structure of the paper is as follows. In Sec. |ll] we present the formulae for the computation of the electron 
spectrum. In Sec. IIIII we describe the numerical evaluation of the spectrum. Section IIVI contains some discussion 
on the different contributions in the high-energy region of the spectrum, and the approximations we have used. We 
conclude in Sec.|Vl where brief comments regarding the implications of our results are given. Appendix [K\ details the 
conventions we use for the Dirac equation, and the electron and muon wavef unctions. 



II. FORMULAE FOR THE ELECTRON SPECTRUM 



The Fermi interaction that mediates muon decay is given by 

Cf = -2V2Gf [V^.^T^^'l^m] [^elpPL^u^] + h.c, (3) 

where Gp = 1.1663788(7) x IQ-^GeV"^ is the Fermi constant and Pl = (I - 75)/2. This La grangian can be Fierz 
rearranged to charge-retention ordering, 

Cf = 2^/2GF [4>eYPLi^t.] [i>.,lpPL^.^] + h.c, (4) 

which is the form that we will use. Since quantum electrodynamics (QED) interactions do not affect the neutrino 
part of the Lagrangian, it is convenient to partition the phase space and integrate the neutrino portion. In that way, 
we generate an effective fi — e current and the free muon decay rate can be written as 



with 
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where the spinors u{p) in that expression are normalized according to u^{p)u'^{p) = 2m5^^ , qp = {Eq,k) is the 4- 
momentum transferred to the neutrinos, Pe is the electron three-momentum and Ee and are the electron and 
muon energies, respectively. When we consider the bound muon decay case, Eq. ([5]) gets replaced by 
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where we are taking the nucleus as static (we discuss the inclusion of recoil effects in the next Section) , and (fe and 
represent the solutions of the Dirac equation for the electron and the muon, respectively. We incorporate the average 
over the muon spin in the definition of while we do not incorporate the sum over the electron spin in the definition 
of <y9e . For the normalization convention that is implied for the wavefunctions in Eq. ([7|) we refer to Appendix [X] The 
muon energy in Eq. ([7]) is given by — — i?b- When the muonic atom is formed, the muon cascades down almost 
immediately to the ground state, the IS wavefunction should therefore be used for </j^ in Eq. ([7|) (the cascade process 
also depolarizes the muons We take the electron to be massless, since electron mass effects are only relevant for 

Ef. ~ TOe, which is not our region of interest^. Integrating over dq^ in Eq. ([7]) we obtain 



E 



% I d'kJ^J^' {q'g,. - gpg.) , (8) 
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where it is understood that Qp = [E^ — E^,k), and we defined 

JP= I d?re-'^--^^,YPL^t.- (9) 



The condition > determines the limit of integration for |fc| to be |fc| < E^ — Ee- Performing the angular integration 

over r in the currents J'', and the angular integrations over k and pe and summing over electron spins in Eq. ([8|) we 
obtain 
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[{E, - E,)k] 2Im( '^^^ 1 + ( 'ToT^T.'^^' + \Sk.\' ] ), (10) 



where 



Glml 

is the free-muon decay rate. The S functions in Eq. (1101) are defined using the notation (...) = Jp°° . . .r^dr, and the 
two cases refer to odd/even + K, respectively 
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{jK+i{kr) [{k - K)f,G +{k + K + 2)g^F]) ' 



^ _U{jK{kr)if,G~g,F)) 
''''~\{jK{kr)ig^G + UF)) ' 

where G and F (<;„ and /«) are the upper and lower components of the radial muon (electron) wavefunction, respec- 
tively. K is the quantum number appearing in the Dirac equation (see Appendix [X] for the conventions we use and 
the definitions of and For a given value of K in Eq. (ITU|) . k can only take the values ±K and zL{K + 1). The 
sum over K goes from to cx), but K cannot take the value K = in the S^^^ and S^^^. terms, and k can never be 
equal to 0. inlz) is the spherical Bessel function of order n. Equation (jlOp agrees with the expressions presented in 
Refs. 



to U. 7ti 



^ In the low-energy region of the spectrum one should also consider the possibility that the electron remains bound or captured by the 
nucleus. 
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A. Inclusion of recoil effects 



In the previous Section we considered the nucleus to be static. The upcoming conversion experiments plan to use 
an aluminum target (previous conversion experiments used heavier elements), where the atomic mass of aluminum 
{Z = 13, A = 27) is 25133 MeV. The nucleus is, therefore, more than 200 times heavier than the muon (m^ = 
105.6584 Me V) and recoil effects should be negligible for most of the electron spectrum. However, the nuclear-recoil 
energy modifies the endpoint of the electron spectrum E^^^, see Eq. ([5]), which means that recoil effects need to be 
carefully considered in studies of the high-energy part of the spectrum. We will always consider the recoil effects at 
first order in a l/m^v expansion, where uin is the mass of the nucleus N. 

For muon DIO, the nuclear- recoil energy is 

-E'rcc = -P^, (13) 



where the three-momentum of the nucleus is 



PN = -Pe - Pu, ~ , (14) 



with ^ and p^^ the three-momenta of the neutrino and anti-neutrino, respectively. We see that, for a given electron 
energy, the nuclear-recoil energy is not constant but depends on the momenta of the neutrinos. This complicates the 
integration over the neutrino momenta, but in the high-energy end of the spectrum we can approximate 

2™ AT 2m N 2mN 2mN^ 

so that the recoil effects amount to a change in the momentum transfer to the neutrinos. The net effect is to substitute 

El , . 

En-Ee^En-Ee —, 16) 

^ ^ 2mN 

in the upper limit of the integration over k and inside the square brackets in Eq. (jlOp . The endpoint of the electron 
spectrum is given by 

E^ E^ 

E^, = E^--^=m^-Eb--^, (17) 
2niN 2mN 

which is exact up to corrections of order 1/m^. The approximation for the recoil energy in Eq. ()15p is the same as 
used in Ref. [l5j |. The electron spectrum including nuclear-recoil effects is, therefore, given by 
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+2k ( E,, - E, 



El \^ jK.{S-^\* + St\*) 
2mN / ™ 2K +1 



.2 / \^k\ + ^kV^ 



+^ V {2K + lf +1^^--' )^ (18) 

where it is understood that this expression should only be used in the region where Eq. ()15p is a good approximation 
to the nuclear-recoil energy. As will be manifest in the following Sections, recoil effects become negligible before 
Eq. ([T5|) ceases to be a good approximation to the recoil energy. That means the inclusion of recoil effects beyond the 
approximation considered here is unnecessary. 



B. Endpoint expansions 



Equation (IT51) constitutes our final result for the high-energy region of the electron spectrum and it is what we will 
use in our numerical evaluations. Nevertheless, it is still interesting to perform a Taylor expansion of Eq. (jl8l) around 
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TABLE I. Values for the B coefficient of tfie leading-order Taylor expansion in Eq. (If 9|) . for a few elements. We use finite-size 
nuclei, characterized by a two-parameter Fermi distribution (see Eq. (|2f p ). with the values of the parameters of that distribution 
taken from Refs. 



Nucleus 


B(MeV-®) 


A1(Z = i3) 


8.98 X iO"" 


Ti(^ = 22) 


4.94 X iO-^'' 


Cu(Z = 29) 


i.i4 X iO~^^ 


Se{Z = 34) 


i.62 X iO"^^ 


Sh{Z = 5i) 


3.57 X iO~^^ 


Au(^ = 79) 


4.79 X iO~^^ 



the endpoint, to make the behavior of the spectrum manifest. We obtain 

r2 \5 



1 dr 
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(19) 



where = {g-^G), = {f-K,F) , r„ = {g-nF), and it is understood that the electron wavefunctions and in 
Eq. ([T9|) correspond to the energy Ee — E^ — E^/{2mN)- We shovif the values of the B coefficient in Eq. ([T9|. for a 
few elements, in Table ID The corresponding Taylor expansion for the expression without recoil effects in Eq. PU)) is 
given by 



1 dr 



To dE, 



64 



{E^ — Ee 



(20) 



where it is understood that the electron wavefunctions in this equation correspond to the energy E,, — E^. Our Taylor 
expansion agrees with the results in Ref. (15| . 



III. NUMERICAL EVALUATION OF THE SPECTRUM 



We now use Eq. ^8)) to obtain a numerical evaluation of the high-energy region of the electron spectrum. We present 
the results for the case of an aluminum nucleus (Al, Z = 13), which is the target intended to be used in Mu2e and 
COMET 

We consider a nucleus of finite size, characterized by a two-parameter Fermi distribution p(r), given by 

pM=Po K^- (21) 

For the parameters of the Fermi distribution we use the values [Tsj 

ro = 2.84 ±0.05 fm, a = 0.569fm. (22) 

Po in Eq. (I2ip is the normalization factor, which can be expressed in terms of tq and a. For the muon mass, aluminum 
mass and the fine-structure constant we use the values = 105.6584 Me V, TOai = 25133 MeV, a = jffQ^, and 
remember that we take the electron to be massless. We numerically solve the radial Dirac equations for the muon and 
the electron, with the charge distribution in Eq. (j2ip . to obtain the wavefunctions. For the muon energy we obtain 



which gives the endpoint energy 



= - = 105.194MeV, (23) 



E^e ^E^- = 104.973 MeV. (24) 
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FIG. 1. Endpoint region of the electron spectrum for aluminum. The squares correspond to the spectrum with recoil effects, 
Eq. I|18p . For comparison, we show the spectrum neglecting recoil, Eq. (|IOp . as the triangles. The right plot is a zoom for 
-Be > 100 MeV, the solid (dashed) line on this plot corresponds to the Taylor expansion around the endpoint with (without) 
recoil. 




FIG. 2. Detail of the electron spectrum for aluminum very close to the high-energy endpoint with (neglecting) nuclear recoil, 
represented as the solid (dashed) line. 



Electron screening will increase the end-point energy in Eq. (1^^ by about -1-0.001 MeV and similarly shift the overall 
spectrum. That small effect is negligible for our considerations. Recall that the sum over K in Eq. ([T5| goes from 
to 00, we include as many terms in K as necessary in order to get three-digit precision for each point of the spectrum. 
This requires about 30 terms near and fewer terms in the low- and the high-energy parts of the spectrum. 

We present the result of the numerical evaluation of the high-energy region of the electron spectrum in Fig [T] 
The squares in the figure are the spectrum with recoil effects, from Eq. (|18p. For comparison, we also show the 
result obtained by neglecting recoil effects, from Eq. pUj) . as the triangles. The right plot in the figure is a zoom for 
Ee > 100 MeV, the solid and dashed lines on this plot correspond to the Taylor expansions in Eqs. ([T9l) and (|20p . 
respectively. Terms up to K = 4 were included in Fig. [T] Figure [5] presents a detail of the electron spectrum very 
close to the high-energy endpoint in linear scale. We can appreciate in that figure how the spectra with (solid line) 
and without (dashed line) recoil effects tend to zero at the corresponding endpoints (the endpoint without recoil is at 
Eg = E^). To make our results easier to use, we mention that the polynomial 

P(£;e) = as^^^ + ae^^ + ay^^ + a85^ (25) 

with 

as = 8.6434 x 10"", ag = 1.16874 x 10"^^ ay = -1.87828 x 10-^^ ag = 9.16327 x 10"^°, (26) 
the energies expressed in MeV, and 

E'^ 

5 = E^^E,--^, (27) 
2toai 



fits very well the result for the electron spectrum in aluminum normalized to the free decay rate (squares in Fig. [T]) 
for all Ee > 85 MeV (i.e., the difference between Eq. (P5|) and the squares in Fig. [T] is not larger than the uncertainties 
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FIG. 3. Electron spectrum for aluminum. Left plot: linear scale; right plot: logarithmic scale. 



discussed in the next Section). Note that, in order to obtain a better fit for the whole Ef, > 85 MeV region, the value 
of as, in Eq. (psj) . was not constrained to be that of the leading coefficient of the Taylor expansion in Table HI 

For completeness, we also show the spectrum for the full range of electron energies in Fig. [3] as the circles, from 
Eq. ()10|) . Terms up to ii" = 31 were included in this plot. The total decay rate for muon decay in orbit in aluminum 
is obtained by integrating the spectrum in Fig. [31 The result we obtain is 

— / —dE, = 0.9934, (28) 

i Jo "-^e 

in agreement with Ref. 14]. Nuclear- recoil effects are negligible in the total rate. Also note, the integrated ordinary 
muon decay rate is hardly affected by the presence of the Coulomb potential, in accord with the general results in 

Sana. 

Since the Mu2e Collaboration also considers titanium (Ti, Z ~ 22) as a viable target [6|, we give the polynomial, 
P^'^''>{Ee), that fits the result for the electron spectrum in titanium (normalized to the free decay rate) for energies 
Ee > 85 MeV, 

p(-)(i?,) ^ + 4-)5f,,) + 4-)<5[,.) + 4-)^f..), (29) 

with 

4^') = 4.44278 X 10"^^ 4^'^ = 9.06648 x 10"^^ a^^'^ = -4.26245 x 10"^^ 4^'^ = 8.193 x 10"^^ (30) 
the energies expressed in MeV, 

^(T.) = E^-Ee- (31) 

and for titanium E^, = 104.394 MeV. 



IV. DISCUSSION 



As expected, the results in Figure [T] show that nuclear-recoil effects are only important close to the high-energy 
endpoint of the spectrum. The corrections to the approximation of the recoil energy that we have used in Eq. (jl5p 
are of order £'e(-E^ — i?e)/(2mjv) (while Ref. jl5l| seems to wrongly estimate this correction as being smaller ^ 
{Efj, — Eef /(2mjv)). For electron energies around 85 MeV, Eq. is still a good approximation to the recoil energy 
while the effect of recoil in the spectrum is very small. When the corrections to the approximation in Eq. (jl5p become 
order one, the recoil effects on the spectrum are negligible. Therefore we conclude that, as anticipated in Sec. Ill A[ 
inclusion of recoil effects beyond the approximation considered here is unnecessary. 

As already noted in Ref. (isj . the Schrodinger wave function for the muon is not a good approximation near the 
endpoint. In that region, one needs to produce an electron with E,, ~ |pe| ~ m^. This implies that, either the 
muon has \p^\ ~ (i-e., it is at the tail of the wavefunction) or (if the muon has the typical atomic non-relativistic 
momentum, of order the inverse Bohr radius) the electron must interact with the nucleus to get \pe\ m^. Those 
two contributions are of the same order in a, which means that we cannot treat the muon within a non-relativistic 
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FIG. 4. Total rate of decay-in-orbit events, for aluminum, with electron energy larger than x, normalized to the free muon 
decay rate To- 

approximation. There are, thus, some leading contributions where the muon is far ofF-shell (it has ^ ^ m^), 
a fact that also tells us that this is the region where finite-nuclear-size effects will be most important (since the muon 
will be closer to the nucleus). By this argument, we can also understand that only the lowest values of the angular 
momentum in the electron wavefunctions contribute at the endpoint (as we see in Eq. (jl9p ). 

Uncertainties in the modeling of finite nuclear-size effects can induce errors in the electron spectrum and, as we 
discussed in the previous paragraph, those are expected to be most important in the endpoint region. We re-computed 
the spectrum varying the parameters in the Fermi distribution as indicated in Eq. (|22p . and found that the errors 
induced in the spectrum do increase as we approach the endpoint, as expected, but they are never larger than ±2%; 
so, we can safely ignore them. 

Finally we comment that radiative corrections have not been included in Eq. (jlSp . However, they are not expected 
to significantly modify the results presented. 

V. CONCLUSIONS AND EXPERIMENTAL APPLICATIONS 

We have performed a detailed evaluation of the high-energy region of the electron spectrum in muon decay in orbit. 
Our results in Eq. Eq. (1^ and Figure [1] include all the relevant effects to accurately describe the high-energy 
region of the spectrum, and provide the correct background contribution for fi — e conversion search experiments. To 
summarize our findings, we plot in Fig. 2] the rate of decay in orbit events producing an electron with energy higher 
than X, as a function of x (normalized to the free muon decay rate). 

The complete muon DIO electron spectrum presented here provides a check on previous low- and high-energy partial 
calculations [iBl ; as well as an interpolation between them. It properly incorporates recoil and relativistic effects 
in the high-energy endpoint region, which is of crucial importance for future fi ~ e conversion background studies. 
In that regard, its primary utility is twofold. First, when experimental data on muon DIO becomes available, our 
formula in Eq. (j25|) can be compared with it and used to refine the detector's acceptance, efficiency and resolution. 
Second, the expected spectrum can be convoluted with the spectrometer resolution function to obtain a more precise 
estimate of the muon DIO background to CLFV fj, — e conversion in the Eg ~ 103.5 — 105 MeV signal region. 

Currently available estimates by the Mu2e Collaboration @ find for 4 x IQ-'^^ ordinary muon captures in Al (corre- 
sponding to a total of 2.6 x 10^® DIO) a signal of 4 conversion events if i?^Ai = T (/i^ Al — > e~Al) /F (/z~Al — >■ ^^^^Mg) — 
10~^^, with only about 0.2 DIO background events for their detector resolution function. At that level, the discovery 
capability is quite robust. However, if i?^Ai is much smaller, more running will be required to enhance the signal and 
refinements of the background using the spectrum in Eq. (|25p will be critical. 

One can also use our analysis to make a rough comparison of DIO backgrounds for stopping targets with different 
Z. Indeed, if a signal for fi e conversion is found in Al, other targets will be important for extracting the underlying 
"new physics" responsible for it. Of course, for much higher Z, the initial dead time of 700 ns envisioned in the Mu2e 
proposal for eliminating prompt backgrounds would significantly reduce the number of "live" muon captures and 
severely compromise the experimental sensitivity. Ignoring that issue for now, we can ask: What DIO background is 
expected for a higher-Z target with energy resolution identical to the Al setup of Mu2e, while continuing to require that 
it also produces 4 signal events? Considering the case of a titanium target, with Z = 22, we expect R^n — l-6i?^Ai 
for models of "new physics" 0] dominated by chiral changing CLFV. So, one needs a "live" run with "only" 2.5 x 10^® 
ordinary muon captures or correspondingly 0.43 x 10^^ total DIO events to reach the same 4 event discovery sensitivity 
as in aluminum. However, even though the total number of DIO events (for all Ee) is smaller by a factor of 6 in Ti, 
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the relative branching fraction for high-energy DIO events in the signal region (with similar detector resolution) is 
about 6 times larger for Ti compared to Al. So, overall the DIO backgroimd is about the same for Ti. 

More difficult will be the loss of muon events in higher Z materials due to the 700 ns dead time during which most 
of the muons undergo capture. For that, a complete reassessment of the muon production and stopping conditions 
may be required. 
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Appendix A: Conventions 



In this Appendix we explain our conventions for the Dirac equation, and the electron and muon wavefunctions. 
The Dirac equation in a central field is given by 



-n5<7r ( ^ + - - ) + ^kr) + TO/3 

or T T 



(Al) 



where W is the energy of the particle, to its mass and Vif) is the potential. The 4x4 matrices 75, /3, K and Ur are 
given by 



75 



1 

1 



,/3 = 



1 
-1 



a ■ r 
B -r 



K 



a-l + l 
+ 



(A2) 



with I the orbital angular momentum I = —if x V and a the 2x2 Pauli matrices 



1 

1 



The wavefunctions are generically denoted as follows 



-i 

1 



1 
-1 



(A3) 



(A4) 



they diagonalize the operators K, p and jz (j being the total angular momentum) with eigenvalues —k, j{j + 1) and 
/U, respectively. Qk, and are the radial functions which are given by the equations 



dfn ^ K-1 
dr r 
dg, 



f.-iW-m-V)g^, 



{W^V + m)U-'^g^. 
dr r 

(f ) are the spin-angular functions, which satisfy 

{a ■ 1 + l)x^ = -KX^ , i.x!: = MX^ / dQfX^Jx''^' = S^^^S, 

They are given by 



(A5) 
(A6) 

(A7) 
(AS) 



with C{lsj;lzSzjz) the Clebsch-Gordan coefficients, Y^^ the spherical harmonics and x™ the spin 1/2 eigenfunctions 

^V = 7X'",SzX™ = mx'",m = ±i, (A9) 
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where s = ^a. Eq. (IA8|) makes manifest that is an eigenfunction of (T • Z + 1 = j ^ — / ^ — s ^ + 1 with eigenvalue 



Thus we have 



I for .7 = ^ — 5 
-/ - 1 for j = ? + i 



(AlO) 



(All) 



and we see that n can take all integer values except 0. We also note that the value of j is given by k according to 



and that the value of I is also given by k, according to 



We express the 15 muon wave function as 



(A12) 



(A13) 



(A14) 



where as is the amplitude of the muon state with spin projection s, since we need an unpolarized muon we have 



\as\ = 1/2. The muon wavefunction is normalized according to 

J r^{F^ + G^)dr = 1. 
We express the electron wavefunction as an expansion in partial waves, according to 

= ^^^^^^^ = E ( -f ^; I ' 

where t is the z-component of the electron spin, and the a^Mt coefficients are given by 



An 



V2 



Uf^^t ^ i'"'—^C [lK-jK;fJ--tt^]Yi'^^ *{pe)e ' 



(A15) 



(A16) 



(A17) 



with Sk the Coulomb phase shift (the distortion from a plane wave due to the potential of the nucleus) . The electron 
wavefunctions are normalized in the energy scale, according to 

(A18) 



where tp'^ ^ corresponds to a solution with energy W. 
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